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Abstract 


The photodissociation dynamics of HBr adsorbed on a LiF(OOl) surface is stud- 
ied using time-independent quantum mechanics. The photodissociation lineshape 
and the Br( 2 P 1 / 2 )/Br( 2 P 3 / 2 ) yield ratio are computed and compared with the cor- 
responding quantities for gas phase photodissociation. The angular distribution 
of the hydrogen photofragments following excitation of adsorbed HBr is computed 
and found to agree qualitatively with experimental data [Bourdon and coworkers, 
J.Chem.Phys. 95, 1361 (1991)]. The effect of polarization of the photon is illustrated 
and discussed. We find the field polarization to affect significantly the magnitude of 
the photodissociation signal but not the angular dependence of the photofragment 
distribution, in agreement with experiment and in accord with expectations for a 
strongly aligned adsorbed phase. 
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1. Introduction 


The photochemistry of surface adsorbed molecules contains potentially a wealth 
of valuable information regarding the orientation, the alignment, and the existence 
of ordering of the adsorbates, the substrate- adsorbate interaction potential, and in 
some cases also the surface structure. 1,2 The study of photodissociation reactions of 
adsorbates is thus of importance for understanding and hence optimizing processes 
such as laser induced surface processing, heterogeneous catalysis and film growth. 3 ’ 4 

Whereas photodissociation of molecules adsorbed on metal surfaces is frequently 
hindered and complicated by competing processes such as quenching and substrate 
mediated excitation, 5 photodissociation of molecules adsorbed on insulators is con- 
trollable experimentally 6-13 and more amenable to theoretical studies. 13- " 0 

Polanyi and coworkers, in a series of time-of- flight mass spectrometry experi- 
ments, investigated the photodissociation of several species of small molecules ad- 
sorbed on clear crystal LiF surfaces. 6-9 Further experiments were carried out by 
Tabares et al., 10 who probed the photodissociation of CH3Br/LiF(100), and by 
Trentelman and coworkers, 11-13 who studied the photodissociation of methyl iodide 
adsorbed on MgO(lOO). The experimental work has stimulated several recent theo- 
retical studies. 13-20 Most of the theoretical work employed classical simulations 13-18 
although a few quantum mechanical models have been also reported. 19,20 McCarthy 
and Gerber, using molecular dynamics simulations, studied the photodissociation of 
ICl/MgO(001) 14 and more recently the photodissociation of CH3l/MgO(001). 13 Strong 
correlation was found between the orientational preference of the adsorbed phase 
and the angular distribution of the photofragments. 13,14 Classical simulations of 
adsorbate photodissociation were carried out also by Huang and Guo, 15,16 by Wat- 
son and coworkers, 1 ' and by Barclay et al. is Guo and Schatz 19 and McCarthy et 
al. 20 reported reduced dimensionality quantum mechanical studies of CH 3 l 19b and 
of IBr 20 adsorbed on a MgO surface. Both groups investigated the effect of the sur- 
face on the electronic branching ratio of the atomic photofragments. The classical 
framework 13-18 has the obvious advantage that motion of the surface atoms can be 
taken into account, at least within a finite cluster model. Classical simulations, how- 
ever, treat the photodissociation event as an instantaneous switch from the ground to 
an excited electronic surface which, in case more than a single excited electronic state 
is involved, may be followed by a curve crossing prior to dissociation. This sequen- 
tial picture, although convenient, is misleading. Quantum mechanical calculations, 
on the other hand, are able to describe accurately the photodissociation dynam- 
ics, but cannot take into explicit account motion of the surface atoms (except via 
much simplified models 21 ) and hence temperature effects are mostly neglected. 19,20 
From these previous studies 6-20 it may be concluded that the presence of the sur- 
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face can affect the photodissociation dynamics in a variety of qualitatively different, - 
system dependent ways, and that adsorbate photodissociation experiments contain 
considerable interesting information that is yet to be explored. 

In the present work we study the TJV photodissociation of HBr adsorbed on 
a LiF(OOl) surface using time-independent quantum mechanics. Several desirable 
features of HBr/LiF have led to our choice of system. LiF has served as substrate 
in most of the experimental photochemistry on insulators studies to date 1,6-10 due 
to its high transparency to TJV light and to its chemical stability. In particular, 
studies 9 of HBr/LiF( 00 l) have shown that the heavy Br atom remains adsorbed 
to the surface subsequent to the dissociation while the hydrogen photofragment is 
ejected into the vacuum. Thus, the difficulty of dealing with a double continuum is 
avoided. Furthermore, the experiments 9 were able to distinguish clearly and char- 
acterize separately hydrogen photofragments that suffered inelastic collisions -with 
the surface or with coadsorbed molecules and those that suffered no such collisions. 
The photodissociation lineshape and angular distribution of the latter fragments 
is of prime interest for the present study since they contain direct information re- 
garding the surface effect on the dissociation dynamics, which is not affected by 
subsequent interactions. Several interesting features of HBr/LiF were noted in pre- 
vious studies : 9,22,23 The angular distribution of the energetic H atoms, the ones 
that dissociated with no loss of energy, was found to peak at a: 55° to the surface 
normal, being independent of the azimuthal (in the surface plane) angle, while the 
Br( 2 ) /Br( 2 P 3 / 2 ) yield ratio was found to be the same as that following gas phase 
photodissociation . 9 Due to the small radius of hydrogen and to the size match of the 
LiF lattice constant and the equilibrium HBr bond length, HBr/LiF(001) exhibits 
rather strong electrostatic attraction; Fourier transform IR experiments 22 and clas- 
sical simulations 23 of the adsorbed phase have shown that the Br atom lies over a 
Li + ion, with the hydrogen atom lying approximately over a neighboring F - ion, 
slightly tilted with respect to the surface plane (1H 0 -116° from the surface normal). 
The photodissociation of HBr/LiF is thus expected to reflect the strong electrostatic 
attraction and orientational preference of the adsorbed phase. 

Classical trajectory studies of the HBr/LiF(00l) photodissociation were carried 
out by Huang and Guo , 15 who reported calculations of the H photofragment angular 
distribution, and more recently by Barclay et al. 18 who emphasized the role of surface 
photodissociation as a flrst stage in surface- aligned chemistry . 1,8,9 Nevertheless, in 
view of the light mass of the hydrogen fragment and the nonadiabatic interaction in 
the excited HBr manifold, a quantum mechanical study appears warranted. 

The goal of the present work is to complement the work of Huang and Guo 15 
and Barclay and coworkers 18 by carrying out quantum mechanical calculations of 
the HBr/LiF( 001 ) photodissociation dynamics. Although the calculations reported 
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below agree with experimental measurements, the major aim of this study is not 
so much to reproduce the insightful observations of Bourdon et a/., 9 but rather to 
emphasize additional aspects and possibly point out directions for future research. 
Our computational method is an extension of a recently proposed time-independent 
scheme for the calculation of gas phase photodissociation dynamics. 24 The method 
was shown in Ref. 24 to be particularly efficient when only the total photodissoci- 
ation cross section is of interest. For practical applications of adsorbate photodis- 
sociation, for instance modeling of chemical vapor deposition reactors, 4 this is more 
often than not the desired quantity. Partial cross sections are required in order to 
compute the angular distribution and branching ratio of the photofragments, and 
their calculation is emphasized in the present study. 

The next section describes our computational model and details the approxima- 
tions involved. In Sec. 2a we describe our choice of potential energy surfaces while 
in Sec 2b we review briefly our technique of computing photodissociation cross sec- 
tions and suggest a method for calculating the photofragment angular distribution. 
Section 3 presents and discusses our results; we compute the photodissociation line- 
shape and the Br( 2 P 1 / 2 )/Br( 2 P 3 / 2 ) yield ratio for gas phase HBr and compare the 
results with the corresponding quantities for surface adsorbed HBr. We compute and 
compare with experiment 9 the angular distribution of the hydrogen photofragments 
and discuss our results in terms of the substrate-adsorbate interaction. Finally, we 
illustrate the role played by the polarization of the light field. Section 4 briefly 
concludes the paper and points out directions for future theoretical work. 
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2. Method 


The application of quantum mechanics to the study of the complex dynamics 
involved in the photodissociation of an adsorbate requires invoking approximations. 
The present model considers the photodissociation of a single molecule, ignoring 
inelastic collisions of the photofragments with the surface and motion of the surface 
ions. The first approximation may be justified by the fact that even in the presence 
of inelastic collisions, their effect can be eliminated from the experimental data, 9 
and the second by the observation that surface relaxation plays a relatively mi nor 
role in the HBr/LiF photodissociation, 15 at least at the experimentally interesting 
temperatures (typically 85 K). Since the heavy Br atom travels on dissociation % 
79 times slower than the hydrogen atom, and was found experimentally to remain 
adsorbed to the surface after dissociation, 9 its motion is neglected in the present 
study. In order to reduce further the dimensionality of the problem the hydrogen 
atom is restricted to move in a plane perpendicular to the surface. It will be the 
goal of future work to systematically relax the aforementioned assumptions and to 
examine their validity more quantitatively. 

Figure 1 shows schematically the configuration of the initial bound state, based 
on the IR measurements of Ref. 22 , and illustrates our choice of coordinate system. 
We use spherical polar coordinates with the origin at the Br atom and the polar 
axis defined to lie along the normal to the surface. In terms of the coordinates of 
Fig. 1 the system (field free) Hamiltonian is given as 


1 1 d r2 9 L_ !L • JL\ 

2 fi R 2 dR OR 2 /iR 2 \ sin 7 $7 ^ d~f ) 


+ V(R, 7) 


(la) 


where /z is the hydrogen mass and atomic units ( h= 1) have been used. With the 
standard choice '$(R, 7) = x{Ri~/)/R the Schrodinger equation takes the form 
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1 d 2 


2 /j. dR 2 
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dx 2 


2x 


d_ 
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-V{R,x)\x(R,x) = 0 (lb) 


a. Potential Energy- 

Having neglected motion of the surface ions, we find that the potential energy 
in Eqs. (1) consists of two terms; the HBr potential energy curves V a [R) and the 
substrate-adsorbate interaction V 3 _ a (R, 7), 

V{R n ) = V a (R) + V,_ a (R, 7). ( 2 ) 

Electronic structure studies of gas phase HBr were carried out by Chapman and 
coworkers, 25 who employed relativistic ab-initio theory to compute the ground and 
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several of the excited HBr electronic states as well as the transition dipole func- 
tions coupling the ground state with the excited manifold. The lowest electronic 
absorption band of HBr extends from « 150 nm to « 220 nm and peaks around 180 
nm. 26 In principle, absorption to three excited electronic states is possible at this 
energy; the 3 Hi and 1 Hi, which correlate asymptotically with ground state Br atoms 
[Br( 2 P 3 / 2 )] &nd are coupled with the ground HBr state via a perpendicular transi- 
tion, and the 3 n o + state which correlates with excited Br atoms [Br* =Br( 2 P 1 / 2 )] 
and is coupled with the ground state via a parallel transition. Chapman ef a/,, 
confirming earlier calculations by Mulliken, 27 showed however, that the transition 
dipole function coupling the ground with the 1 III state is over an order of mag- 
nitude larger than those coupling the ground state with the 3 13i and 3 H 0 + states. 
The experimental detection of ^ 14% excited Br atoms 28 has been attributed to 
spin-orbit coupling between the initially excited 1 IIi and the 3 n o +^ 25 A fourth- ex- 
cited electronic state, the 3 Ei, which correlates with Br( 2 P 3 / 2 ), is accessible at the 
relevant energies through nonadiabatic interaction. Nevertheless, since the coupling 
with this state becomes significant only at much larger H-Br separations 25 its effect 
on the dissociation dynamics is expected to be small. In the present work we thus 
restrict attention to three (diabatic) electronic states; the ground ( 1 E 0 -f) state, the 
1 IIi state which is coupled with the 1 E 0 + by the electromagnetic field, and the 3 n 0 + 
state which is coupled with the 1 Hi by spin-orbit interaction. The relevant potential 
energy curves, following a fit to the data of Ref. 25b, are shown in Fig. 2. We fit 
the dipole moment function of Chapman et al. to the form 

f i(R) = D d ex V [-b d (R-R <i ) 2 ] (3) 

with D d = 0.38 a.u., b d = 1.2 Bohr -2 and R d = 1.8 Bohr. A functional form 
for the spin-orbit coupling potential was derived from analysis of the electronic 
wavefunctions of Refs. 25. We used the asymmetric form 


V so (R) = D c 


exp[— b el (R - .Re) 2 ] 

1 + exp [—b e2 (R - Rc)] 


( 4 ) 


with D c = 3.34 x 10 -3 a.u., b cl = 0.3 Bohr -2 , b c2 = 1.5 Bohr -1 and R c = 4.5 Bohr. 

The substrate-adsorbate potential V 3 - a in Eq. (2) consists of two terms: an elec- 
trostatic term, denoted Vf_ a , which accounts for the interaction between the electric 
field generated by the surface ions and the charge distribution of the molecule, and a 
core term, denoted V^_ oJ which describes the short range repulsion of the closed shells 
and the long range attractive dispersion interactions. The electrostatic forces can in 
principle be computed by straightforward summation over the pairwise Coulombic 
interaction of the lattice ions with the adsorbate. This procedure is, however, very 
difficult to apply in practice due to the long range of the Coulomb force (in fact the 
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sum of Coulombic terms is known to be only conditionally convergent). As pointed 
out by Lennard-Jones and Dent 29 and by Steele, 30 the electrostatic potential $ gen- 
erated by a regular lattice of surface ions can be accurately expressed as a rapidly 
converging Fourier series. For the (001) face of an FCC crystal 30 


+9l )/ 2 


4e _ , 

= . ; 

“■ a * \J3i - 9z | 


exp 


27 rz‘ 


- \fgl + dl 


x cos 


27 r 


X 

9 1- 


y 

92 - - 

a 


1 + exp 
9i + 92 


-*yi 

)] 


9l 


92 


]} 


( 5 ) 


where e is the electron charge, a is the lattice constant (a = 4.03 A for LiF) and 
the prime signifies summation over odd values of g 2 and g 2 . x and y in Eq (5) are 
measured in the surface plane, along the [100] and the [010] crystallographic axes, 
respectively, and z' denotes distance along the normal to the surface, measured with 
respect to the positive ion chosen as the origin ( z ' = z + i2u- Br in Fig. 1). For many 
physisorbed systems the adsorbate equilibrium position is sufficiently far from the 
surface ( z ' > a) that only the leading terms in Eq. (5) (|< 7 i| = |^| = 1) contribute 
significantly to the sum. Since, as mentioned in the introduction, the HBr molecule 
lies relatively close to the surface, more terms should be included but, as is clear 
from the exponential z' dependence of $(r), the sum converges rapidly. Following 
Ref. 23 we model the molecule as a distribution of three collinear point charges 
with qi = 0.276e located at the H atom, q 2 = — 0.444e located at the Br atom and 
q z = 0.168e located 1.292 A from the Br end of the molecule, on the opposite side to 
the H atom. This choice is not unique; 31 other models, such as a two point dipoles 
model, 23 are possible and were investigated [vide infra). With the above choice the 
electrostatic substrate-adsorbate interaction is written as 


VU = £ »*( r,) 

k-\ 


( 6 ) 


and thus depends on the molecular orientation relative to the lattice axes. Since 
the excited HBr electronic states are repulsive and the molecule dissociates into 
neutral atoms the electrostatic interaction is omitted from the excited potential 
energy surfaces. 

The core interaction V^ c _ a , assumed pairwise additive, can be readily computed 
by direct summation over the surface ions. Following Refs. 18 and 23 we use a 
Tang-Toennies form 

V TT {q) 

f 2n(<?) 




fc=o 


kl 


( 7 ) 
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where the first (Born-Mayer) term accounts for the short range repulsion and the 
second term is a damped dispersion series which is typically truncated after its 
second term. The various parameters in Eq. (7) are taken from the semi-empirical 
calculations of Barclay and coworkers. 18 

The substrate-adsorbate interaction potential is shown in Figs. 3 as a function of 
height above the surface (z') and distance along the [100] direction of the crystal (x). 
Figure 3a corresponds to the bound state, where corrugation is more pronounced and 
the electrostatic part of the interaction introduces dependence on the orientation of 
the molecule. Figure 3b corresponds to atomic (free) hydrogen, where the interaction 
is periodic in x and considerably less corrugated. 

Using Eqs. (2, 4-7) the complete HBr/LiF(001) potential energy is given as 

V g (R, 7 ) = u 12o ji?)+y;_ a (i7,7) + K c -a(^7) 


V e (J2, 7) 


( K Ul (R) + v;i{(R, 7 ) 
V v 30 (r) 


v> 


3 n, 


0 + 


V, 0 (R) 

(R) + v;i{(R,i) 


( 8 ) 


where g stands for ground, e stands for excited, and V7il a and denote, respec- 
tively, the core interaction [Eq. (7)] in the bound and in the free states. 


b. Photodissociation Dynamics 

With the field free Hamiltonian fully specified, the photodissociation cross sec- 
tion may be computed via a similar method to that described in Ref. 24. The total 
cross section 


= H^n/4( w ) 


n ,1 


'4tt 2 


U) 


was written as 


cr;(u>) = - 


/4 irw\ 


) 


\ C 

47ru>\ 


lM-e| *i,o)| 

n>l 


Im(«i i0 \p ■ £ G{E + ) p ■ e\ 


^ Im ^2 ( Xi /x 0 ,/Gi,/-(^ + )^,o|xi) 


( 9 ) 


( 10 ) 


where $ ii0 ) is a bound state of the matter Hamiltonian, is an incoming 

scattering state of the same Hamiltonian, the indices i and n denote collectively the 
set of rovibrational quantum numbers, and l is an electronic index, p • e in Eq. (9) 
is the component of the transition dipole moment in the direction of the field, the 
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ground state wavefunction in Eq. (10) is written as a product of an electronic state 
and a nuclear state, ^i >0 } = | 0 ) | Xi), Ho ,l = ( 0 | p. ■ e | /), and Gi t i>(E + ) is the 
matrix element of the outgoing Green’s operator 

G(E + ) =]jm(E + ie- H)- 1 (11) 

in the electronic basis. 

While the total photodissociation cross section of Eqs. (9, 10) reflects the short 
time dynamics, partial cross sections are required in order to compute observables 
that depend on the long time dynamics, such as the branching ratios and quantum 
state or angular distributions of the photofragments. Using the Lippmann-Schwinger 
equation for the scattering states in Eq. (9), 32,33 

n»j) = [i + £(£-)!'''] | 

= — , ( 12 ) 

where | ^ are solutions of an approximate Hamiltonian Hi = H — V', the 

partial cross section can be expressed formally as 24 

= (^T-j \eG(E + )fl-k\ $i, 0 }| 2 

= H xi) ■ (is) 

We used absorbing boundary conditions (ABC) 24,34 to obtain a well behaved rep- 
resentation of the Green’s operator; thus, the infinitesimal convergence factor e in 
Eq. (11) was replaced by a coordinate dependent operator e(q) which is practically 
zero in the interaction region but becomes sizable outside the physically relevant 
region of space. A discrete variable representation (DVR) 35 was employed for the 
nuclear degrees of freedom. 

In application of Eqs. (12, 13) to surface photodissociation, the boundary con- 
ditions that are to be satisfied by the scattering solutions are clearly different from 
those appropriate to gas phase photodissociation. We first note that the method 
takes into exact account energy conserving collisions of the photofragment with the 
surface (i.e., collisions that do not transfer energy to the surface, resulting in the two 
high energy peaks in the energy distribution of Ref. 9). This can be seen by expand- 
ing the Green’s operator in terms of its gas phase analogue G a (E+) = (E+ie — H a )~ l 
where H a = T +V a [see Eq. (2)]. Substituting the expansion in Eq. (10), for example, 

= - (^) Im {( *» I G.(£ + ) | <S„ ) 

+ (*,! G.(£ + )V)_.G.(E+) !*„) + ...} (14) 
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where we denoted | ^ ) = /2-e | 

_ir S|( $ M I $ ^,;)| 2 > where 

n,Z 


^io y The first term in Eq. (14), Im ( | G a (E + ) | 4^ ) 

|^£?n,z) are eigenstates of H ay corresponds to the 


direct photodissociation. Subsequent terms incorporate the effect of multiple, en- 
ergy conserving collisions with the surface. This is an essential feature for the 
present system since, as pointed out in the introduction (and discussed in previous 
studies 9,22,23 ), the HBr molecule is well oriented in the adsorbed phase, with the H 
atom tilted toward the surface and hence redirecting collisions of the photofragments 
with the surface are expected to play a role. 

Below we will be particularly interested in calculation of the photofragments’ 
angular distribution, i.e., the differential photodissociation cross section <r(k, u). 
To that end we first write the Lippmann-Schwinger equation in the coordinate rep- 
resentation as 


( 


R I 


E,k,l 



+ JdR'G(E--,R,R’)V(R') (r' | 


(15) 


where k is a unit vector in the direction of observation and we indicate explicitly the 
dependence of the scattering wavefunction on both the magnitude and the direction 
of k. The right hand side of Eq. (15) vanishes in the crystal ( z — » — oo) and takes 
the standard asymptotic form of a scattering solution 32 at large separation from the 
surface ( R — r oo). In Eq. (15) and below we omit the n index from the scattering 
waves and confine attention to the case of a single displacement vector (here the 
H-Br displacement vector R) of interest for the present application. The method 
is, however, equally applicable to the dissociation of polyatomic adsorbates as will 
be evident from the discussion below. A simple choice for the zero order solutions 
are (energy normalized) plane waves, 



where ki = — V^°°) and V t 00 is the asymptotic potential in the Zth electronic 

state. Using the resolution of the plane wave 32 


e' kR = 4* 


( 17 ) 




where jj(kR) are spherical Bessel functions, 36 we obtain 

R I 


E,k,l 


l2fik[ 


7T 


m(k) r) 


J,TTl 
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I <LR' G(E ~ ; R, R') V(R') 


(18) 
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where V’ of Eq. (15) has been replaced by the complete interaction V. With the 
notation 

{ R | fl, m ) = iAkR)Yi.^m (W) 

the scattering state is written as 

) = E YUV [i + g{e-)v\ | ) (20) 

j,m 


and, using Eq ( 12 ), the differential photodissociation cross section assumes the form 

i 1 2 


<^,i(k,U>) 



E^w(a°„^ 

j,™ 


eG{E + )p. -k 



( 21 ) 


We note that 
V = 0, thus 


/°jm) satisfy Schrodinger’s equation in spherical coordinates for 

J dkcr^k,^) = <T lt i(u>) 


where is given in Eq. (13) and we used the orthonormality of the spherical 

harmonics, fdkY? m (k)Yji i7n >(k) = SjjsSm^. For systems where the initial state 
is aligned at an angle 7 ^ < 7t/2 with respect to the normal, as is the case for 
the photodissociation of CHsI , 11 ^ 13 for instance, Eq. ( 21 ) predicts that the angular 
distribution of the photofragments will peak at an angle 9 k & 7 The symmetry 
of the spherical harmonics in Eq. (21) ensures that in systems such as HBr/LiF, 
where 7 ; > tt/2 and the photofragments are likely to collide with the surface after 
dissociation, a(9k^) will peak at an angle close to that expected for specular recoil 
( 0 jt ^ 7 T — 7 {). Deviation from the specular direction which depends on the degree 
of surface corrugation, and a certain width which reflects the degree of alignment, 
are expected. 

In practice one would normally replace the plane wave solutions by distorted 
waves, 


'I' 1 ~ 
EXi 




( 22 ) 


where the potential has been partitioned as V = V\ + V' and G\(E~) is the Green’s 
operator (modified to include an absorbing potential) corresponding to H\ = T + Vi. 
Thus 



[i + G(£-)V"]|* 1 a,) 

- E **„£) <?(£-) < <?,(£-)« I ) 


-iY,YUi)G(E-)t 



(23) 
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where are defined by the last equation. Equation (12) 

our choice of zero order Hamiltonian, thus | ) can be any 

able) basis. Below we choose 


holds regardless of 
(physically accept- 


( R I fl.j.m ) = 9(R COS 7) ( R I ) 


(24) 


where g(z = R cos 7 ) is a cutoff function which vanishes as the wavefunction pen- 
etrates into the crystal and approaches 1 elsewhere. For more complicated appli- 
cations it would be desirable to employ a better zero order basis, for instance by 
including a part of the physical gas-surface interaction in the distorted Hamiltonian, 
so as to reduce the computational effort. 
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3. Results and Discussion 

The photodissociation lineshape of gas phase HBr, computed via Eq. (10), is 
shown as a dashed curve in Fig. 4. The overall structure and the peak position 
agree with the experimental data of Lee and Suto. 26 The percentage yield of excited 
Br atoms [Br* =Br( 2 P 1 / 2 )], obtained by employing Eq. (13) to compute the partial 
Br and Br* photodissociation cross sections, is « 15% of the total Br yield. The 
solid line in Fig. 4 shows the photodissociation lineshape for surface adsorbed HBr. 
The lineshapes have the same overall structure but the one corresponding to the 
adsorbed phase is slightly distorted and blue shifted by 2.2 nm. We found the ground 
state to be stabilized in the presence of the surface by ~ 0.18 eV as compared to 
the gas phase. (A somewhat smaller energy shift, ~ 0.15 eV, was found with a 
two dipole model for the electrostatic interaction 23 with a corresponding slightly 
smaller shift of the photodissociation lineshape.) The stabilization, a significant 
part of which is due to electrostatic interaction, is related to the adsorption energy 
of HBr. At a given photon energy the translational energy of the photofragment 
is thus lower by ~ 0.18 eV from the H-Br relative translational energy following 
gas phase photodissociation [2.6 eV for the Br( 2 P 3 / 2 ) channel and 2.14 eV for the 
Br( 2 P x / 2 ) channel]. The decreased translational energy of the H photofragment 
was found also in the experiments of Bourdon and coworkers 9 and in the classical 
trajectory simulations of Barclay et a/., 18 although the reverse effect was found in the 
classical work of Ref. 15. Since the photodissociation cross section is proportional 
to the photon frequency (in the weak field limit) the energy shift of the lineshape is 
somewhat smaller than the stabilization energy of the ground electronic state. The 
Br*/Br yield ratio was found to be unaltered from its gas phase value, as was found 
also experimentally. 9 This result, however, is to be expected since the form of the 
coupling potential and the relative position of the excited electronic states in the 
present model are taken to be the same as in the gas phase. 

In Fig. 5 we show the angular distribution of the energetic [Br( 2 Pz/2) channel] 
hydrogen photofragments following excitation of HBr/LiF(001) at 193 nm. The 
distribution of H fragments corresponding to the Br* channel was found to have the 
same structure but a smaller overall magnitude. As shown in Fig. 5, the distribution 
peaks at ~ 56° with respect to the surface normal, a smaller angle than that expected 
for specular recoil, i.e., for scattering from a flat surface (69° — + 64° for an initial 
orientation of 111° — ► 116° with respect to the normal). The effect of surface 
corrugation is thus to tilt the scattered photofragments toward the surface normal. 
Figure 6 shows the experimental angular distribution of Bourdon and coworkers 9 for 
three different surface coverages. The peak position of the measured curve (~ 55°) 
and its structure at large angles are reproduced well by the calculated curve of Fig. 5, 
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although the latter falls off more rapidly at small angles. This discrepancy may be 
expected since the model considers a single adsorbed molecule while the experiments 
were carried out at a finite coverage (0.2 — > 1.0 monolayers). Figure 6 shows clearly 
that the effect of increased coverage is to increase substantially the small angle wing 
of the distribution while leaving the large angle contribution practically unaltered. 
As discussed in Sec. 2 the present calculation, which neglects motion of the surface 
ions, reflects most realistically the limit of zero temperature, while the experiments 
of Bourdon et al were carried out at 85 K. It appears, however, that our neglect 
of surface relaxation is not the reason for the discrepancy between the calculated 
and the experimental distributions at small angles since the present distribution 
is similar (although broader) to the calculations of Ref. 18. The latter, based on 
classical simulations and a single excited electronic state model, considered likewise 
an isolated adsorbate but allowed for temperature effects. 

The angular distribution of the fragments is expected to depend on the adsorbate 
orientation, the effect of redirecting collisions and the polarization of the field. The 
adsorbate orientation determines both the direction at which photofragments will hit 
the surface and the orientation of the transition dipole function with respect to the 
space fixed frame. The field polarization affects the dissociation probability through 
the dot product /x-e in Eqs. (9, 13, 21). The electronic transition under consideration 
is perpendicular 26 and hence jl'k — /x(i?)e*£ where e is a unit vector perpendicular to 
R and /x(i?) is given in Eq. (3). The effect of the initial orientation of the adsorbed 
phase is contained in the form of the ground state wavefunction, which, in turn, 
reflects the substrate-adsorbate interaction in the ground electronic state. Simple 
Franck-Condon type considerations based on Eq. (21) indicate that in the case where 
the initial state is strongly aligned, 0 has a large amplitude for a narrow range 
of 7 values which give rise to dominance of certain partial waves ({j,m}). The 
latter determine the structure of the angular distribution. In that situation the field 
polarization is expected to affect strongly the magnitude of the signal but have little 
effect on its angular dependence. If, on the other hand, the initial state samples a 
wide range of orientations, different polarizations will accentuate the contribution of 
different scattering components and hence alter the form of the angular distribution. 
Similar conclusions were reached in Ref. 9 based on different considerations. In 
Fig. 7 we compare the computed angular distributions following dissociation with 
light polarized along the x and the z axes. As expected, the signal corresponding 
to x polarized light is considerably smaller than that corresponding to z polarized 
light (light polarized along the normal to the surface) while the overall structure of 
the angular distribution is practically invariant to the polarization. Figure 7 is in 
agreement with polarized light measurements of Bourdon and coworkers 9 who found 
the angular distribution of the H fragments to be nearly independent of the field 


14 



polarization. This result provides an independent support of the finding that the 
HBr molecule is strongly aligned in the adsorbed phase. A crude estimate of the 
preferred orientation angle may be obtained by comparing the magnitudes of the 
x and the z polarized light signals. We find cr(k,a>;e = x)/cr(k,w;£ = z) » .38 
indicating a preferred orientation angle of 7 l 111 0 with respect to the normal, in 

qualitative agreement with experimental IR 22 and photodissociation 9 results, and 
with molecular dynamics simulations of the adsorbed phase 23 ( 7 ; ~ 111 0 — 116°). 

It should be remarked that the structure of the angular distribution depends 
quite sensitively on the substrate-adsorbate interaction potential, as may be ex- 
pected in view of the above discussion. Nevertheless, the agreement of our calcu- 
lated distribution with experiment cannot serve to verify the quantitative details 
of the potential function of Eqs. (5-7) but only its qualitative features, namely, its 
corrugated nature and its manner of orienting the initial state. 
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4. Conclusions 


In this work we developed a time-independent, quantum mechanical model for 
study of adsorbate photodissociation and applied it to the photodissociation of HBr 
adsorbed on a LiF(OOl) surface. The photodissociation lineshape of adsorbed HBr 
was found to be blue shifted with respect to that of the gas phase molecule, consistent 
with stabilization of the ground electronic state in the presence of the surface. The 
computed yield ratio of excited Br atoms [Br* =Br( 2 Pi/ 2 )] does not change from its 
gas phase value according to the present model, as was found also experimentally. 9 
The theoretical angular distribution of the H photofragments peaks at ~ 56° to the 
surface normal in agreement with experiment (55 db 5 0 ) 9 although the experimental 
distribution extends to smaller angles. In view of the trend of the experimental 
distribution with increasing coverage, it is likely that the discrepancy is due to -the 
presence of coadsorbates in the experiment which is not taken into account in our 
calculations. The role played by the held polarization was examined and discussed. 
We found the polarization of the electromagnetic field to change significantly the 
magnitude of the signal but not its angular dependence, in agreement with polarized 
light experiments 9 and in accord with qualitative expectations for a strongly aligned 
initial state. 

As pointed out in the introduction, existing theoretical studies of adsorbate pho- 
todissociation can be broadly classified as either purely classical 13 ” 18 or simplified 
quantum mechanical models. 19,20 The former are able to take into account surface 
relaxation but treat the photodissociation event and the effect of nonadiabatic cou- 
pling in the excited manifold in an approximate manner. The latter, to which the 
present model belongs, are able to treat accurately the photodissociation process 
but are forced to neglect motion of the surface ions (or else introduce crude approx- 
imations for the interaction). The theory of adsorbate photodissociation therefore 
appears ideally suited for application of system-bath type approximations 3 ' where 
the system degrees of freedom are treated by exact quantum mechanics and the bath 
modes are approximated. It allows, moreover, a natural sepaxation of the degrees of 
freedom into an active subspace and a remaining, passive subspace. In particular, 
like other chemical processes at surfaces, photodissociation reactions of adsorbates 
have aspects in common with liquid phase processes, where the solvent is treated 
statistically . 38 Doren and Tully 3 ' a studied the role played by precursor states in ad- 
sorption and desorption dynamics using a mean field approach to cast the effect of 
all modes but the reaction coordinate mode in the form of a temperature-dependent, 
effective potential for motion along the reaction coordinate. Of the wide variety of 
system-bath type approximations which have been developed in recent years 37 the 
approach of Ref. 37a appears most attractive since it introduces a time-independent 
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effective potential for motion of the system, and since a large number of bath modes 
can be taken into account. Ideally one would treat statistically only the surface 
modes, and when applicable, the modes of coadsorbates. Nevertheless variations 
that include in the bath several of the adsorbate modes are likely to be needed (and 
valid) in photodissociation studies of larger molecules. Refinement and extension 
of the present model via this and other approaches will be the objective of future 
research. 
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Figure Captions 


1. Schematic illustration of HBr adsorbed on a LiF surface following the IR mea- 
surements of Ref. 22. 

2. Potential energy curves for HBr (derived from Ref. 25). 

3. Substrate-adsorbate interaction as a function of the Cartesian coordinates of 
the H atom, z' denotes height above the surface and x denotes distance along 
the [100] crystallographic axis with Li + defining the origin. The Br atom is 
fixed at its equilibrium position. 

(a) for bound (molecular) HBr 

(b) for free (atomic) H+Br. 

Contours are in increments of leV. 

4. Total photodissociation cross section for HBr as a function of the photon 
wavelength: 


( -) gas phase HBr 

( ) surface adsorbed HBr. 


5. Angular distribution cr(9k, A = 193 nm) of the hydrogen photofragments cor- 
responding to the Br( 2 P 3 / 2 ) channel following excitation of HBr/LiF(001) at 
193 nm. 9k is measured with respect to the surface normal. 

6. Experimental angular distributions of the hydrogen photofragment correspond- 
ing to the Br( 2 P 3 / 2 ) channel at different surface coverages: 

( o o o o ) 0.2 monolayer 
( • • • • ) 0.5 monolayer 
( q o o □ ) 1.0 monolayer. 

(Reproduced from Ref. 9 with permission of the authors.) 

7. Angular distribution of the H photofragments following excitation at 193 nm: 

( ) light polarized along the z axis 

(- ) light polarized along the x axis. 
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